DNA in nanopore - counterion condensation and coion depletion 
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Abstract 

Molecular dynamics simulations are used to study the equilibrium distribution of monovalent 
ions in a nanopore connecting two water reservoirs separated by a membrane, both for the empty 
pore and that with a single stranded DNA molecule inside. In the presence of DNA, the counterions 
condense on the stretched macromolecule effectively neutralizing it, and nearly complete depletion 
of colons from the pore is observed. The implications of our results for experiments on DNA 
translocation through a-hemolysin nanopores are discussed. 

PACS numbers: 87.15.Aa, 87.15.Ya, 87.14.Gg 
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Recent experiments on translocation of a single strandedfss) DNA molecule through an 



a-hemolysin nanopore inserted in a membrane 



theoretical studies 
through narrow pores [8|, 




5| have prompted a number of 



. These studies, as well as other models of polymer translocation 
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14 lisj l , focused on the polymer aspects of the 



problem such as the effect of confinement-induced entropy losses on the translocation of 
a Gaussian chain and on the effective friction between the chain and the narrow pore. 
The fact that DNA is a charged object was considered only as far as the mechanism of 
pulling by an externally applied potential gradient was concerned. While the above theories 
were quite successful in reproducing the observed distribution^f translocation times ^ and 
the molecular weight dependence of the translocation velocityjTi , some important questions 
remain unanswered. For example, as has been noted in reference 6|, the theory overestimates 
(by more than an order of magnitude) the effective pulling force on the DNA. A similar 
discrepancy concerning the magnitude of the force acting on DNA has been observed in 
recent experiments on nanopore-induced opening of DNA hairpins 0. Furthermore, while 
the models focus on the kinetics of DNA translocation through the pore, experiments probe 
this process only indirectly, by monitoring the transient blocking of the current of small 
ions through the channel. It is therefore essential to understand the physics of these ions 
in confined space - their interactions with the pore and with the DNA molecule, as well as 
their mutual interactions. Notice that even though the inner diameter of the a-hemolysin 
pore is much larger than the size of a small ion, the presence of negatively charged DNA is 
expected to have a dramatic effect on the distribution of small ions in the pore. 

We proceed to study the problem by molecular dynamics simulations. We take an isolated 
rectangular box of three sides 56 A, 56 A and 112 A in the x, y and z directions, respectively. 
This box is separated into upper and lower compartments by a membrane of thickness 50 A. 
The center of the membrane is permeated by a cylindrical pore of radius 7.5 A positioned 
along the central z-axis and assumed to be neutral (we do not account for the presence of 



charged and polar amino-acids at the membrane- water interface jl2|). In the absence of DNA 
the simulation box is filled with three particle species, namely, counterions, colons and water 
molecules that are free to move throughout the upper and lower compartments and through 
the pore, but can not permeate the volume occupied by the membrane. The counterions 
and colons are spheres of radii 1.3 A and 1.8 A corresponding to ionic radii of and Cl~ 
respectively, having either a positive or negative unit charge at their centers. There are 100 



2 



salt ions of each sign in the simulation box, corresponding to salt concentration of about IM. 
Water molecules are emulated by putting a sphere of radius 1.2 A in every (2.9 Ay volume, 
typically amounting to 6800 particles (a value smaller than the Van der Waals radius of water 
was taken in order to speed up the simulation). The electrostatic properties of water and 
of the membrane have been taken into account by treating them as dielectric continua with 
Eu, — 80 and £^ = 2 respectively, and matching the two dielectrics by introducing a 2.1 A 
thick transition layer at the water-membrane boundary, across which the dielectric constant 
changes linearly from the former to the latter value. Reflecting boundary conditions are 
assumed at the surfaces of the rectangular box and of the membrane, and at the pore wall. 

The ssDNA is modeled as a flexible chain of 12 repeat units each of which contains a 
backbone made of a spherical bead with charge — e at its center (phosphate group) connected 
to a non-charged bead (sugar ring), to which an additional non-charged bead (base) is 
attached. The beads are connected by effective springs and undergo thermal fluctuations 
under the constraint that the center of the macromolecule is pinned down at the pore center. 
The radii of the charged, non-charged, and side-chain beads are taken as 2.1 ^4, 1.8 ^4 and 
2.2 A, respectively (the latter value corresponds to an adenine base). When DNA is present, 
12 additional counterions are added to the system. 

For each instantaneous spatial configuration of the charges, the electrostatic potential 
$(x) obeys the Poisson equation V • (eV$) = — 47rp, where the charge density p is ob- 
tained by summation over ions p(x) = J^- qiS{x — x^), with S being a smeared 5— function 
(with the boundary condition $ = on the surface of the box). This equation is solved 
numerically using the conjugate-gradient method with grid spacings 0.7 A in the x,y- 
directions and 1.1 A in the z-direction. The force on the i-th particle is then given by 
Fj= — V(gi$+f/Lj) + A Qj^ij / ^i^ijj where the last term accounts for the short-range 

correction to the electrostatic force (the summation is taken over nearby ions within a Bjer- 
rum length, 7 A, and is computed in every time step). Since the rapid variation of the 
electrostatic force between neighboring charges is taken into account by this term, the re- 
maining contribution to the electrostatic potential ($) is calculated by solving the Poisson 
equation only once in every 10 time steps. The volume exclusion of particle cores is incor- 
porated through the Lennard- Jones potential t/Lj(rj) = X^^- 4eLj((cr/rij)^^ — (cr/rjj)®) for 

^ 2^/^(7 and Ulj — —^lj otherwise, where — ri — Vj and a is the sum of the radii of 
the two interacting particles. The temperature T is used as the energy scale and €lj — ksT 
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is assumed for simplicity. 

We first consider a nanopore without DNA. We run the simulation for 4000 ps (our time 
step corresponds to 3.4 fs) and average the results over time. Inspection of Table 1 shows 
that on the average there are about 2.2 counterions and the same number of colons in the 
pore; however, the fluctuations are of the same magnitude as the average. This value is lower 
than that (4) corresponding to a uniform distribution of ions in the system, in agreement 
with the expectation that ions are repelled from the pore because the electrostatic self- 
energy of an ion in the pore is higher than in the bulk, due to the low dielectric constant 
of the surrounding membrane jisl. In order to clarify the underlying physics, in Fig. la we 
plot time histograms of the numbers of counterions (A^+), colons (A^-) and the total charge 
{Q) in the pore. Notice that even though the fluctuations around the mean values are 
quite large, these numbers fluctuate nearly exactly in phase and therefore the total charge 
exhibits only small variations (±0.7e) around zero. One is tempted to interpret the above 
observation as evidence for formation of neutral Bjerrum (coion/counterion) pairs^^ inside 
the pore and indeed, such pairs and even longer string-like aggregates of alternating colons 
and counterions) are clearly observable in snapshots of the instantaneous distribution of ions 
in the pore. 

The time histograms of a pore with DNA inside are presented in Fig. lb. The time 
averaged results for the numbers of counterions, colons and of DNA charges in the pore 
are summarized in Table 1. In the physical 2, inspection of Table 1 and Fig. 

2a reveals that while the counterion concentration in the pore reaches nearly 3 times its 
value in the bulk, colons are nearly completely depleted from it. The number of counterions 
appears to be determined by the requirement that the negative charges of the DNA are 
fully neutralized. Similar counterion enrichment and somewhat smaller colon depletion are 
observed for abasic DNA (with bases removed) which has a smaller excluded volume but 
the same charge as normal ssDNA. 

In order to elucidate the role of membrane electrostatics on the distribution of ions in the 
pore, we increased the dielectric constant of the membrane to its value in water, Sm = 80. 
While a 4-fold increase in the concentration of colons in the pore compared to the — 2 case 
is observed, the number of counterions is nearly unchanged (see Table 1 and Fig. 2b). We 
conclude that while the effect of increased electrostatic self-energy of charges surrounded 
by low dielectric constant medium plays a major role in colon depletion, the number of 
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counterions in the pore is determined mainly by the condition of electroneutrahty. This 
leaves open the question of whether the counterions are condensed (i.e., strongly localized) 
on the DNA or free to move around the pore. While it is difficult to distinguish between 
condensed and uncondensed counterions in a narrow pore, inspection of Figs. 2a and 2b 
suggests that counterions are more strongly localized on the DNA charges in the low Em 
case. In order to make a quantitative comparison we calculated the appropriate radial 
distribution functions and found that the probability to find a counterion within 4 A from 
a DNA charge is nearly 2 times higher in the = 2 than in the Sm = 80 case. This is 
consistent with the notion that counterion condensation on DNA is more pronounced in the 
former case, because of the larger electrostatic self-energy of DNA in the pore 18]. Since 
stronger DNA-counterion coupling should lead to stiffening of the macromolecule as the 
results of its "dressing" by the counterions, one expects the persistence length of DNA to 
increase with decreasing Sm- Indeed, we find that the end-to-end distance is larger and that 
transverse fiuctuations are smaller when = 2 than when Sm = 80 (44 A versus 42 A, 
respectively). Further evidence for the renormalization of persistence length scenario comes 
from comparing the end-to-end distances of abasic, uncharged (with charges removed from 
spheres representing the phosphate groups) and normal ssDNA. As bases are added to the 
charged backbone, R increases from 35 to 44 A and when charges are removed from the 
latter R decreases from 44 to 40 A, in accord with the expectation that the addition of 
the bulky bases has a larger effect on the persistence length of DNA than its "dressing" 
by the small counterions. Since nearly complete colon depletion from the pore is observed 
only when both (a) DNA is present and (b) Sm = 2, it must be due to a combination 
of excluded volume and electrostatic effects. A tentative explanation is that while in the 
empty pore case Cl~ ions can enter the pore by forming Bjerrum pairs with K"^ ions, such 
bulky ion pairs are effectively excluded from the pore in the presence of DNA. Since a similar 
mechanism should apply even if the pore is occupied by a neutral macromolecule, we ran the 
simulation with uncharged DNA. We find that there are about 1.4±1.0 colons (and the same 
number of counterions) in the pore, a number intermediate between that corresponding to 
the empty channel (2.2) and the pore with a charged macromolecule (0.4) cases. This again 
is consistent with the expectation that the charged DNA is "dressed" by the neutralizing 
counterions which increase its excluded volume. 

In summary, when ssDNA is present in the pore, it is neutralized by counterion conden- 



5 



sation; coions, on the other hand, are expelled from the channel. These conclusions remain 
valid (A^+ = 13, = 1.5, N^^^ = 12, R = 39 A) even if we increase the pore radius to 10 
A, a value that is closer to the mean radius of an a-hemolysin channel. In spite of the various 
assumptions inherent in the present simulation (treating the solution and the membrane as 
dielectric continua, neglecting hydration, using a simplified model for DNA and taking no 
account of the detailed molecular structure of the pore wall), the emerging physical picture 
of nearly complete counterion condensation and colon depletion in the pore appears to be 
quite reasonable. Analogy with electrophoresis suggests that the condensed counterions fa) 
reduce the effective charge of DNA and (b) move together with it in an electric field^, 0|. 
If this analogy is correct one is tempted to conclude that these counterions do not take part 
in charge transport across the pore and that the experimentally observed drop of the ion 
current in the presence of DNaQ] should be attributed to the reduction in the number of 
charge carriers due to DNA-induced depletion of coions. However, since the identity of the 
char ge c arriers is sensitive to molecular details of the pore neglected in our work (see, e.g., 
ref. [17]), the tentative conclusion that residual conduction through the DNA-blocked pore 
is due to the negatively charged coions should be tested experimentally e.g., by changing 
the chemical composition of the salt molecules. 
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Table 1: The average numbers (and variances) of counterions (A^+), colons (A^-) and of 
DNA charges (N^^"^) in a pore of 7.5A radius embedded in membranes of low (sm — 2) and 
high = 80) dielectric constants, are summarized. The net charge in the pore in units of 
electron charge {Q/e) and the end-to-end distance (i?) of ssDNA within the pore are also 
shown. Only ions and DNA charges with centers further than 1.8 ^4 away from the ends of 
the pore were counted in order to minimize end effects. 
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FIG. 1: (a) The time history of the numbers of counterions (A^+), coions (N^) and the total 
charges {Q = e{N-^ — N-)) contained in the empty pore (no DNA) is shown, (b) The time history 
in the presence of DNA (with n£^^ 12) in the pore. The total instantaneous charge in the pore 
is Q = e{N+ -N-- N^^^). In both cases Sm = 2. 
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FIG. 2: Snapshot of the ion distribution with ssDNA in the pore for = 2 (a) and £m = 80 (b). 
Green and red spheres are counterions and coions, respectively and orange spheres represent the 
charged phosphate groups of DNA (the neutral sugars and bases are not shown) . 
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